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Abstract 

In the paper, several concrete examples, as well as their numerical simulations, are 
given to show that parameter identification based on the so-called adaptive synchronization 
techniques might be failed if those functions with parameters pending for identification in 
coupled systems are designed to be mutually linearly dependent or approximately linearly 
dependent on the orbit in the synchronization manifold. This failure might be emergent 
not only when the synchronized orbit is selected to be some sort of equilibrium or some 
sort of periodic oscillation, but also when it is taken as some type of chaotic attractor 
produced by driving system. This result implies that chaotic property of driving signal is 
not necessary to achievement of parameter identification. The mechanism inducing such a 
failure, as well as the bounded property of all trajectories generated by coupled systems, 
is theoretically expatiated. New synchronization techniques are proposed to rigorously 
realize the complete synchronization and parameter identification in a class of systems 
where the nonlinearity is not globally Lipschitz. In addition, parameter identification are 
discussed for systems with time delay. 



PACS: 05.45.Gg, 05.40.-a, 87.10.+e. 
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1 Introduction 



The most classical phenomenon in reference to synchronization is owing to Huygens' obser- 
vation about the synchrony of pendulum clocks [Ij. Since this historical discovery, synchro- 
nization as an omnipresent technical issue has become a focal topic of great importance in 
many applications. Moreover, the basic concept related to chaos synchronization in cou- 
pled chaotic systems was initially introduced by Pecora and Carrol in 1990 [2]- Since their 
seminal paper, chaos synchronization as an interesting research topic of great potential 
application has been widely investigated and consequently applied in plenty of fields, rang- 
ing from secure communications to pattern recognitions, from complex network dynamics 
to optimization of nonlinear systems, and even from chemical reaction to brain activity 
analysis In particular, a wide varieties of synchronization approaches, including tradi- 
tional linear or nonlinear feedback coupling, impulse coupling, invariant manifold method, 
adaptive design coupling techniques, and white-noise-based coupling have been fruitfully 
proposed [l]-[5] and several types of synchronization, including complete synchronization, 
generalized synchronization, phase synchronization, and lag synchronization, have been 
introduced in succession [6]- [9]. 

Among all the proposed coupling approaches for realization of complete synchroniza- 
tion between coupled chaotic systems with or without time delays, the newly developed 
adaptive design coupling technique has aroused a great amount of attention from many 
researchers |10j-|18j simply due to the reported success in unknown parameter identifi- 
cation. Their explorations have shown that unknown parameters could be identified in 
several well-known chaotic systems and even in some neural network models with or with- 
out time delays. In particular, consider an n-dimensional nonlinear system described by 

x = F{x,p), (1) 

where a; = {xi,X2, . . . ,Xn)^ G K", F{x,p) = (^Fi{x,p), F2{x,p), . . . , Fn{x,p)^ , and 

m 

Fi{x,p) = Ci{x) + ^Pijfij{x), i = l,2, ...,n. (2) 

Here, Ci{x) and fij{x) are, respectively, assumed to be some kind of real valued functions, 
and p = {pij} U C M" are (n • m) parameters pending for identification, in which U 
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is some bounded set. Given the bounded driving signal x{t) produced by system (1), its 
response system is designed through 




Qij 



Sijeifijiy) 



(3) 



z = 1,2,... ,n, j 



= 1,2,... ,m 



where the feedback couphng e • e is in the form of (eiei, £262, . . . , ene-n)^ , &i = — Xi), 
q = {qij}, and both and 5ij are arbitrarily chosen positive constants. 

A question naturally arises: "Is it possible to accurately identify all the {n ■ m) pa- 
rameters of the chaotic system provided that the output time series of system ([T]) are 
experimentally obtained?" When the response system is designed as ([3]), the answer to 
this question, as mentioned above, is reportedly positive. Concretely, the complete syn- 
chronization between the driving system ([T]) and the response system ([3]) could be always 
achieved; moreover, the varying parameters g in ([3]), initiating from arbitrary values, will 
be asymptotically convergent to the correct values of the parameters p in ([1]) as time tends 
towards positive infinity. Seemingly, their theoretical arguments are based on a delicate 
design for Lyapunov function, on the well-known Lyapunov Stability theorem, and even 
on the LaSalle invariance principle. 

As a matter of fact, those signals produced by these driving-and-response systems 
could be completely synchronized; nevertheless, the parameter identification might be 
failed if those terms with parameters pending for identification are designed to be mutually 
linearly dependent or approximately linearly dependent on the synchronized orbit. In the 
paper, not only concrete examples with their numerical simulations will be provided to 
illustrate such a failure of parameter identification, but also the mechanism inducing this 
failure will be anatomized. Also, the performed analysis will show that chaotic property of 
synchronized orbit in synchronization manifold ( or say, chaotic property of positive limit 
set of driving signal) is not always necessary to achievement of parameter identification. 

The rest of the paper is organized as follows. In Section 2, three concrete examples, 
as well as their numerical simulations, are consecutively given to illustrate the possible 
occurrence of parameter identification failure. The synchronized orbits in those examples 
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are, respectively, selected to be some sort of equilibrium, periodic oscillation, and chaotic 
attractor. The mechanism inducing this failure, as well as the the bounded property of all 
trajectories generated by the coupled systems ([T]) and Q, is theoretically expatiated in 
Section 3. Furthermore, in Section 4, new synchronization techniques are further proposed 
to realize complete synchronization and parameter identification in a class of polynomial 
systems where the nonlinearity is not globally Lipschitz. In Section 5, parameter identi- 
fication are discussed for systems with time delay. Finally, the paper is closed with some 
concluded remarks. 

2 Examples Showing Failure of Parameter Identification 

In this section, three groups of driving-and-response systems are concretely presented to 
illustrate the possible occurrence of failed parameter identification. 
First, consider the Lorenz system: 

Xl = Pi{x2-xi), 

X2 = P2X1 - X1X3 - X2, (4) 

X3 = X1X2 - P3X3 
as a driving system. And the corresponding response system becomes: 

yi = qi{y2-yi) + ei{yi-xi), 

m = 92^1 -yiys- 2/2 + £2(2/2 -2:2), (5) 
m = yiy2- qsys + esivs- X3), 

where the updating laws of q = {qi,q2,Q3) and e = (ei,e2,e3) are, respectively, taken 
as: qi = -6i{yi - xi){y2 - yi), q2 = -^2(2/2 - 2:2)2/2, 93 = -hiya - X3){-y3), and 
ei = -n(2/i - xi)"^, €2 = -^2(2/2 - X2f, 63 = -r3(y3 - X3f. 

Particularly, when the parameters are taken as pi = 35, p2 = ^, and p3 = 28, the com- 
plete synchronization between systems (jj]) and ([5]) could be achieved with time evolution, 
which is numerically shown by Fig. 1(a). This group of parameters, which are different 
from the classical parameters inducing strange attractor of the Lorenz system, simply 
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make the synchronized orbit generated by system ^ become an asymptoticaUy stable 
equihbrium E, as is shown by Fig. 1(b). If the reported analytical results are completely 
correct, it could be expected that the varying parameters q will be eventually convergent 
to the accurate values of the parameters p = ipi,P2,P3)- However, this is not the case. 
As depicted in Fig. 2, although the numerical simulation is consistent with the expectation 
for the parameters q2,3, it is beyond the expectation for the parameters qi. Concretely, 
the value of qi does not approach but always keeps a distant from the accurate value 
of pi with time evolution, which means failed parameter identification does occur for qi. 
Intuitively, there must exist some mechanism inducing such differences between qi and 
Q2,3 when adaptive synchronization techniques are taken into account. 

Secondly, construct a driving system based on the Chen's system through: 

Xi = 30{X2 - Xi) + Q{xi,X2,Pl,P2), 

X2 = (28 - 30)xi -xiX3 + 28x2, (6) 

±3 = XiX2 - 3X3. 

where the additional term 

(xi cos(0.9026) + X2 sin(0.9026)) ^ 



Q{xi,X2,Pl,P2) = 0.1x<^pi- 



-P2 



23.442 

, 2 



xi sin(0.9026) + X2 cos(0.9026)) 

7I92 ^ 



and both pi and p2 are parameters expected to be identified. As a matter of fact, without 
the term Q, system ([6|) becomes the original Chen's system admitting an attractive peri- 
odic orbit. As displayed in Fig. 3, the projection of this attractive periodic orbit into the 
X1-X2 plane is approximately looked upon as an ellipse. Thus, when pi = 1 and p2 = —1, 
the term Q actually is an approximate formula of this projection in the X1-X2 plane. 

Given a driving signal (xi, X2, xs)"^ generated by system ([6]), the corresponding response 
system is designed to be in the form of 

yi = 30(y2 - yi) + Q(yi, 2/2, 91,92) + ei(yi - 2:1), 

m = (28 - 30)yi-yiy3 + 287/2 + £2(^2 -X2), (7) 
2/3 = 2/12/2 - 31/3 + £3(^3 - X3). 



in which, according to ([3]), the updating laws of the two varying parameters are taken as: 

{yi cos(0.9026) + 7/2 sin(0.9026)) ^' 
23.442 ' 

( - yi sin(0.9026) + 2/2 cos(0.9026)) ^ 

and the adaptive techniques of couphng strengths are, respectively, chosen as: ei = 
-ri(yi - xi)2, £2 = -r2(y2 - 2^2)^, and 63 = -r^iys - xa)^. 

Contrary to the expectation, both qi and q2, starting from almost every initial values, 
fail to approach the real values of the parameters pi = 1 and p2 = —1, as shown in Fig. 4. 
This example, as well as the first example, shows that failure does occur for parameter 
identification based on the adaptive synchronization technique when the synchronized 
orbit is particularly selected to be some type of steady dynamics, such as asymptotically 
stable equilibrium and attractive periodic orbit. 

Instead of the above-mentioned steady synchronized orbit, the existing numerical re- 
sults |10j-|18j show that parameter identification could be always achieved when those syn- 
chronized orbits in the synchronization manifold are designed to be some type of chaotic 
attractor in advance. Then, a question arises: "Is chaotic property of synchronized orbit 
in synchronization manifold necessary to achievement of parameter identification based 
on the adaptive techniques?" 

To find out an answer to this question, consider a 4-dimensional model developed from 
the original chaotic Lorenz system as a driving system: 

xi = pi{x2-xi), 

±2 = P2X1 - X1X3 - X2, 

(8) 

X3 = XiX2-p3X3+P4:X3{l+xl), 

Xi = 0x4 + — X3), 

where pi = 10, p2 = 28, and = | are three classical parameters for the original 
Lorenz system to generate chaotic attractor, and a = —100, b = 0.1, p4 = 1. Given 
these parameters, the orbit produced by system ([8]) in the synchronization manifold still 
exhibits chaotic property in the phase plane, which is displayed by Fig.5(a)-(b). This 

7 



qi = -Siiyi - xi) 
q2 = -S2{yi-xi) 



chaotic property is further verified by calculating the largest Lyapunov exponent of the 
system (Ai ?a 0.54274 > 0), as is shown by Fig.5(c). 

Provided with the driving signal produced by system ([8]) , the complete synchronization 
between systems (j8]) and its response system could be numerically achieved as long as the 
response system is designed through: 

yi = qi{y2-yi) + ei{yi-xi), 
m = g22/i - 2/12/3 - 2/2 + £2(2/2 - 2:2), 

(9) 

2/3 = 2/12/2 - 932/3 + 942/3(1 + 2/4) + £3(2/3 - 2:3), 

2/4 = a2/4 + ^(2/1 - 2/3) + £4(2/4 - 2:4), 
in which the updating laws of the parameters are taken as qi = —5i{yi — xi){y2 — 2/i)> 92 = 
-^2(2/2 - X2)yi, 93 = -^3(2/3 - X3){-y3), 54 = -^4(2/3 - 2^3) [2/3(1 + , and the adaptive 
coupling strengths are taken as ei = — ri(yi — 2:1)^, 62 = — ^^2(2/2 — 2:2)^, £3 = — '"3(2/3 — 2:3)^, 
€4 = —r4{y4 — 2:4)^. In spite of the success in complete synchronization and in parameter 
identification for qi 2, it is impossible to utilize (73^4, initiating from almost every points, 
to identify the accurate values of the parameters ^3^4 in system ([8]) . All these are shown 
in Fig. 6. Clearly, this example implies that the answer to the above-posed question is 
negative. 

Remark 2.1 The fourth-order Runge-Kutta scheme is used to solve all the ordinary 
differential equations in our numerical simulations. 

3 The Mechanism Inducing the Failure 

On the one hand, three concrete examples in the last section show that some parameter 
identification might be failed no matter what kind of dynamical phenomenon is displayed in 
the synchronization manifold. On the other hand, many existing numerical results always 
show successful parameter identification. In order to clarify the mechanism inducing such a 
seeming paradox, we perform a more delicate argument by adopting the LaSalle invariance 
principle [19] as follows. Similar to [11], set a Lyapunov function candidate by 

7X 7% 771 Tl 

Vie, e, = ^ E + ^ E E ^(%- - P^^)" + '2 E ^(^^ + ^)'- (10) 

i=l i=l j = l ^-^ i=l * 
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Then, the derivative of the function V{e, e, q) along with the coupled systems ([T]) and ^ 
could be estimated by 

n 

V{e,e,q) ^ (n/-L)^ef. 

j=i 

Here, it should be pointed out that / is not the locally Lipschitiz constant of the function 
Fi{x,p) but the uniformly Lipschitiz constant since the bounded property of the trajec- 
tory y{t) generated by the newly response system ([3]) are not confirmed but pending for 
confirmation yet. 

Now, we contend that e{t), e{t), and q{t) are bounded for all t ^ to, where to is the 
initial time. Indeed, one of the three variables is supposed to be unbounded on [to, +00), 
so that V{e{t),e{t),q(t)) is also unbounded on [tQ,+oo) according to (|1U|). On the other 
hand, V{e{t), e{t), q{t)) ^ V{e{to), e(to), q{to)) simply due to V{e, e, q) ^ for sufficiently 
large L. This contradiction thus implies the bounded property of e(t), e(t), and q{t) for 
all t ^ to- 

Therefore, in light of the LaSalle invariance principle, the trajectory (e(t), e(t), 
initiating from any location in the phase plane will eventually approach the largest invari- 
ant set A4 contained in the set 

£ = {ie,e,q) \ y(e,e,q) =0}. 

Then, the main concern becomes how to make a clear description of the invariant set M 
ever contained in the set £ with respect to systems ([I]) and ([3]). To this end, a combination 
of systems ([1]) and ([3]) yields 

m m 

ei = Xi - iji = Ci{x) - Ci{y) + '^Pijfij{x) - '^qijfijiy) - ei{yi - Xi) 

i=i i=i 

m m 

= Ci{x) - Ci{y) + ^ [pijfij{x) - qijfij{x)] + ^ [qijfij{x) - qijfij{y)] - ei{yi - Xi). 
j=i j=i 

(11) 

Also, notice that V{e, e,q) =0 implies e = x — y = 0, ei{t) = 0, and qij{t) = 0. It thus 
follows from ([TT]) that, for every orbit (e(t), e(t), G f of the coupled systems, 

m 

K -%(i)]^-(a3(t)) = 0, (12) 
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where each qij{t) is identical to some constant q*j. And the largest invariant set contained 
in £ with respect to systems ([T]) and ([3]) is 

M = {{e,€,q) I e = x-y = 0, = e*, qij = q*^}, 

where x = x{t) is the synchronized orbit, or mathematically say, the positive limit set of 
the driving signal. Thus, the question becomes: "Is each q^- surely equal to Pij?" From 
(jl2p , the answer to this question is theoretically positive provided that [LIM] : for any given 
i, {fij{x),j = 1,2, .. . ,m| are linearly independent on the synchronized orbit x = x(t) in 
the synchronization manifold. 

For an accurate definition of linearly independent or linearly dependent functions, refer 
to [20]. Also, it is valuable to mention that two functions might be linearly independent in 
a domain but linearly dependent in some subset contained in this domain. For example, 
functions gi{s,u) = s and g2{s,u) = li} are obviously linearly independent in but they 
are linearly dependent in a parabola-like subset 5^ = u) G I s = iiu^] C for 
some nonzero constant [i. 

If hypothesis [LIM] is not satisfied, for some i = zq, either there exist two nonzero 
functions fif^j^(x) and fioj2{^) linearly dependent on the orbit x{t) in the synchronization 
manifold, or simply f^-^j^{x{t)) = 0. We focus on the former case since failure of parameter 
identification could be easily illustrated in the latter case. Accordingly, fi^^j^[x{t)) = 
cfi^j^ixit)) for some nonzero constant c, which at most implies that [pioii ~ %oii(^)] + 
^\_Pioj2 ~ ^«oj2(*)] = 0- Clearly, although qioji{t) and qi^^jiit) are, respectively, identical to 
some constants q*^j_^ and q*^j^, it is not certain that q*^j_^ = pi^j^ and g*^^^ = pjoi2- Actually, 
they are totally distinct in most cases. Therefore, parameter identification might always 
be failed if hypothesis [LIM] is not strictly satisfied in the design of driving-and-response 
systems. 

Next, by virtue of the argument performed above, the reason why parameter identifi- 
cation fails in the three examples given in the previous section is expatiated as follows. 

For the driving-and-response systems dH) and ([5]) with specific parameters, the syn- 
chronized orbit x*{t) = {x\{t) , X2{t) , x'^{t))'^ in the synchronization manifold, as shown 
in Fig.l, is a globally asymptotical equilibrium E = x*{t) = (6.8313,6.8313,1.6667)"^. 
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Substitution of dH) into ^2il gives 



\pi - qiit)Mit) - xlit)] =0, [p2 - q2it)]xlit) =0, [pa - g3(t)]x^(t) = 0, 

where each qi{t) is identical to some constant in the invariant set A4 {i = 1,2,3). 
According to [20], each x*{t) (7^ 0) is hnearly independent and X2{t) — xl{t) (= 0) is 
hnearly dependent. This imphes that 3 is identical to p2,3 but ql is not necessarily 
identical to pi. Therefore, qi{t), though obeying the updating law, will not be surely 
convergent to pi. This illustrates the reason why parameter identification succeeds for 
'72,3 but always fails for qi as shown in Fig. 2. However, when the synchronized orbit 
x*(t) with the classical parameters is chaotic, X2(t) — xl{t) is nonzero and thus is linearly 
independent, which satisfies hypothesis [LIM]. Hence, qi{t) will be convergent to pi almost 
surely, as is shown by many existing numerical results. In addition, when the synchronized 
orbit x*{t) unfortunately becomes the unstable equilibrium of the original chaotic systems, 
X2{t) — xl{t) is still identical to zero violating hypothesis [LIM]. So, qi{t) still will not be 
surely convergent to the accurate value of pi in such a case. 

For the coupled systems © and the orbit x*(t) in the synchronization manifold, 
as mentioned above, is designed to be some kind of stable periodic orbit. Its projection 
into the X1-X2 plane, which seems like an ellipse, could be approximately expressed by the 
formula 

(xi cos(0.9026) + X2 sin(0.9026)) ^ _^ ( - xi sin(0.9026) + X2 cos(0.9026)) ^ ^ _ ^ 



23.442 7.192 
Also, substitution of ([6]) into ([12]) yields 



, (xi{t) cos(0.9026) + X2{t) sin(0.9026))^ 
bi - Hit)] ^2 + 

r , , J (-xi(t)sin(0.9026) +X2(t)cos(0.9026))' 1 ^ 
[-P2 + q2{t)] < ^ YJ^2 - U = 0. 

Thus, as long as the complete synchronization between systems ([6]) and ([7|) is achieved, 
the orbit x{t), as well as y{t), will approximately approach the stable periodic orbit x*{t). 
Then, both functions 

(xt(t)cos(0.9026) +4(t)sin(0.9026))^ ( - (f) sin(0.9026) + cos(0.9026)) ^ 

23AA^ ^"""^ 7T9^ ^ 
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are approximately linearly dependent. This, according to the argument performed above, 
means that both qi and q2 are not suitable for parameter identification, as is verified by 
the simulation results in Fig. 4. 

Unlike the steady dynamics exhibiting in synchronization manifold in the previous two 
examples, the synchronized orbit x*{t) generated by the driving system ([S]) is deliberately 
designed to be chaotic in the sense of possessing positive Lyapunov exponent. Analogously, 
substitution of ([8]) into (fT2|) produces 

[P3 - q3{t)]x3{t) + [p4 - qi]x3{t) {1 + [X4{t)f} = 0. 

It is obvious that functions X3 and 2:3(1 + X4) are linearly independent in the whole phase 
plane M^; nevertheless, they are approximately linearly dependent on the synchronized 
orbit x*{t) because the cubic term [x4(t)]^ is almost equal to zero as time t is sufficiently 
large (see Fig. 7). Thus, this illustrates the reason why 53 and 94 initiating from a mass of 
points will not be convergent to the accurate values of and p4, respectively, in concrete 
numerical simulations. 

In addition, consider a case that parameter b in both systems dH) and ([9]) is selected 
to be zero instead of 0.1. This case could be regarded as a very special example where 
parameter identification may fail in spite of existence of chaos. In such an example, 
because of x|(t) = 0, functions X3 and 2:3(1 + 2:4) are definitely linearly dependent on 
the corresponding orbit x*{t), which violates hypothesis [LIM]. Therefore, (73 and (74 can 
not be utilized to identify the parameters and p4. In a word, chaotic property of 
synchronized orbit in synchronization manifold does not always guarantee a success in 
parameter identification. 

Remark 3.1 In the last two examples, those functions on the orbits in the synchroniza- 
tion manifold are approximately linearly dependent. Rigorously, they are still linearly 
independent in a mathematical sense, so that parameter identification might be theoreti- 
cally achieved for qi correspondingly with pi. However, in real application, discretization 
techniques, such as the Runge-Kutta scheme and the Euler scheme, are always taken into 
account in solving the coupled continuous differential systems. Thus, owing to the preci- 
sion limit, it is unavoidable that dynamics produced by the discretized system may not 
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be completely consistent with the true dynamics generated by the original system. It is 
the approximate dependence of those functions that poses some trap of local critical point 
for Qi and that leads to a failure of parameter identification in the last two examples. 
Therefore, not only rigorous linear- dependence of those functions with parameters pending 
for identification on the synchronized orbit but also approximate linear- dependence on the 
synchronized orbit should be always avoided whenever adaptive synchronization techniques 
are used in practical parameter estimation and chaos communication. 

4 Complete Synchronization without Globally Lipschitz Con- 
dition 

In the previous section, it is pointed out that hypothesis [LIM] is indispensable for a suc- 
cessful parameter identification. In addition, the uniform Lipschitz condition for F{x,p) 
is also important in the argument performed above for obtaining a non-positive property 
of V{e,e,q). As a matter of fact, this uniform condition could be loosed if the bounded 
property of the response system ([3]) could be priorly estimated. However, this prior esti- 
mation could not directly follow from the bounded property of driving system ([1]) since 
dynamics of its response system with coupling term might be completely different from 
the driving system which is, though, supposed to be bounded in advance. Then it poses a 
question: "Other than the above coupling technique and uniform Lipschitz condition, un- 
der what kind of coupling methods and conditions on F[x,p) can one obtain a successful 
parameter identification rigorously?" 

To this end, it is first assumed [HPT]: each Fi{x,p) is homogeneous polynomial with 
degree no more than two with respect to x. As a matter of fact, a large quantities of 
nonlinear systems does not satisfy globally Lipschitz condition but are consistent with 
this assumption [HPT], such as the Lorenz system and the Chen's system. 

Next, notice that 

2yfcyj - 2xkXj = {yk - Xk){yj + xj) + {yj - Xj){yk + Xk) 

— Gj I j G I 2^ G-j 

13 



for arbitrary k and j. Then, it is easy to verify that each ei[Fi{y,p) — Fi{x,p)~^ can be 
written as a homogeneous polynomial with degree no more than three with respect to 
e = y — X if assumption [HPT] holds. 

Reasonably, the driving signal x{t) generated by system ([T]) is supposed to be bounded 
in advance. In order to obtain a rigorous synchronization in the system where the nonlin- 
earity only satisfies assumption [HPT], we re-designed the response system as: 

y = F{y,q) + e-e + u -e^, 

ii = -Tief, uii = -Sief, (13) 

Qij ~ ^ij^ifij{y)i 



where u ■ = (wief , u;2e2, • • • jW^e^) , each Si is arbitrarily positive constant, and other 
states and parameters are the same as those defined in ([3]). 
Set a Lyapunov function candidate by 




Thus, the derivative of this function along with the coupled systems ([TJ and yields 

n n n 

Hie, e, u, q){t) = ^ ei{t) [Fi{y{t),p) - Fi{x{t),p)\ - ^ Mel{t) - ^ iVef (i), 

j=l i=l i=l 

where both M and N are positive numbers. By virtue of the conclusion on each Cj [Fi{y, p) — 
Fi{x,p)'j obtained above, the elementary inequality 

l=i,j,k 

and the assumed bounded property of the driving signal x{t) and parameter set U, we can 
obtain that H{e,e,u:,q){t) ^ for sufficiently large numbers M and N. 

By using the similar argument performed in the previous section, we can easily prove 
that every trajectory generated by the coupled systems ^ and ([13]) is not only bounded 
for all t ^ tQ but also approaching the largest invariant set contained in 

£" = {{e,e,u:,q) \ H{e,e,u,q) = O} 
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with respect to these coupled systems. More precisely, the largest invariant set becomes 

M' = {(e,e,cj,q) | e = 0, = e*, tOi = uj*, qij = g* }, 

where e*, lo*, and q*j are some constants dependent on the initial values of the coupled 
systems. Furthermore, to achieve an accurate parameter identification between systems ([T]) 
and (fT3]l . hypothesis [LIM] should be still adopted. Then, the above performed argument 
could be concluded as the following proposition. 

Proposition 4.1 If assumptions [LIM] and [HPT] on F{x,p) are satisfied, the complete 
synchronization between driving system UP and its response system il'j]) could be surely 
achieved, and the parameter identification could be accurately realized in a mathematical 
sense. 

Remark 4.2 As mentioned above, in numerical experiment and even in real application, 
not only hypothesis [LIM] should be strictly satisfied but also the approximate linear- 
dependence attributed to precision limit should be avoided. 

Remark 4.3 Assumption [HPT] on F{x,p) could be further generalized to some other 
case where globally Lipschitz condition is not fulfilled. For instance, one could further 
consider the system where the degree of the homogeneous polynomials is larger than two, 
or some of those polynomials are non- homogeneous. However, additional coupling terms 
(e.g. e^^^^ , f = 2, 3, • • • ) should be added into the response systems in order to obtain a 
successful synchronization and parameter identification in a rigorous sense. 

Remark 4.4 As a matter of fact, nonlinearities in the previous three examples are not 
globally Lipschitz but polynomial. Theories and coupling techniques (i.e. Proposition 4.1 
and Remark 4.3) proposed in this section should be utilized to deal with those systems 
for obtaining a successful synchronization and parameter identification. 

5 Parameter Identification in Systems with Time-Delay 

Time-delay, as an omnipresent phenomenon, can not be neglected in practice. So, in this 
section, complete synchronization and parameter identification in time-delayed systems 
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via adaptive coupling techniques are further investigated. For simphcity, consider a one- 
dimensional driving system: 

x{t)=af{x{t))+bg{x{t-T)), (14) 

where r ^ is a time-delay, a and b are parameters pending for identification, and 
functions / and g are assumed to be globally Lipschitz continuous with Lipschitz constants 
kf and kg, respectively. Given driving signal x(t) generated by system (dH), the response 
system is designed as 

y{t) = a{t)f{y{t)) + p{t)g{y{t - r)) + r,{t)e{t) + uj{t)e{t - 5), 

ait) = -f{y{t))e{t), (j{t) = -g{y{t - T))e{t), (15) 

fj{t) = -e^(t), u;(t) = -e{t)e{t - 6), 

where (5 ^ is a time-delay induced by coupling term, error dynamics e{t) = y{t) — x{t). 
The initial conditions for coupled system (jl4p and (jlSp are chosen as e = (p, a = A, P = B , 
T] = E, CO = W £ C = C{[— maxjr, 5}, 0], M), in which C denotes the sets of all continuous 
functions from [— max{r, 6}, 0] to M. 

Set a Lyapunov functional candidate by 

V{cP,AB,E,W) = ^cp\0) + ^[A{0)-af + l[B{0)-bf + \[E{0) + Lf 

where L, M are some proper positive constant. Then, the derivative of V along with 
coupled systems and ([13]) could be estimated by 

V{<P,A,B,E,W) ^ (^akf + 2-^'^cf>''{0)+bkg-\m\-\H-r)\-4>H-r) 



Clearly, V{4', A, B, E, W) becomes non-positive provided L > max <j + 2akf -|- 



M2 b^'kl 



4 J . By using a similar argument performed above, we can conclude that every trajec- 
tory (et(i;^>), Of (A), /3i(-B), ?7((£'), cJt(VF)), starting from arbitrary initial condition, is surely 
bounded for all t ^ — max{r, S\. 
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Then, according to the invariance principle for the systems with time-delay |21j . every 
trajectory, as time tends towards positive infinity, approaches the largest invariant set A4 
contained in 



with respect to coupled systems (|14p and (llSp . This further implies that the first com- 
ponent of each element in Ai is identical to zero (i.e. = 0) and the others are some 
constant functions (i.e. A = A* , B = B*, E = E*, and W = W*). The accurate values of 
these constant functions rest on the initial conditions of coupled systems (I14p and ()15p . 

Parameter identification could be achieved if both equations A* = a and B* = b 
are valid. However, these equations are not always tenable although = indicates a 
successful complete synchronization between systems ()14p and (jlSp . As a matter of fact, 
subtraction of (jl4p from the first equation in (jl5p yields, in A4, 



which follows from e{t) = y{t) — x{t) = in M. Now, it is clear that when functions 
f{x{t)) and g{x{t — r)), on the synchronized orbit x{t) in the synchronization manifold, 
are linearly dependent, both A* = a and B* = b are not certainly tenable. More concretely, 
(i) when the driving signal asymptotically tends towards some equilibrium x{t) = x* of 



linearly dependent so that parameter identification for a and b will be almost surely failed; 
(ii) when the synchronized orbit x{t) is periodic with period r and both functions / and g 
are linearly dependent in R, parameter identification also will be failed; (iii) when x{t) is 
chaotic, parameter identification will be achieved in a mathematical sense for non-constant 
differential functions / and g, and even for f = g (see an example shown in Fig. 8(a) 
where both / and g are taken as sinusoid functions). However, for case (iii), parameter 
identification also might be failed in numerical simulation or in real application. For 
example, in spite of chaotic property of x{t), it is likely that x{t) ~ x{t — r) with a small 




= a{t)f{y{t))-af{x{t))+m9{y{t-T))-bg{x{t-T)) 



[a{t)-a]f{x{t)) + [l3{t)-b]g{x{t-T 



)) 



[A*-a]f{x{t)) + [B*-b]g{xit-T)) 



system ([HI, two constant 
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time-delay or that fluctuation of x{t) seems to be relatively steady in a macro scale. These 
extraordinary cases may lead to approximate linear-dependence between functions /(^(t)) 
and g[x{t — r)), which thus results in failure of parameter identification in numerical 
simulation. See an illustrative example in Fig. 8 (b). In addition, r = could be regarded 
as a special case where parameter identification is always failed provided that functions / 
and g are linearly dependent on x{t). 

In conclusion, we have the following proposition on synchronization and parameter 
identification for coupled systems p4|) and psp with time-delay. 

Proposition 5.1 The complete synchronization between driving system and its re- 
sponse system U5\) could be surely achieved via adaptive coupling techniques. Furthermore, 
the parameter identification could be accurately realized in a mathematical sense provided 
that f{x(t)) and g{x{t — t)) are linearly independent on the synchronized orbit x{t) in the 
synchronization manifold. 

Remark 5.2 With analogous arguments but more complicated notations, the results on 
the driving system ()14|) could be further generalized to the case where higher dimensional 
driving systems and multiple parameters identification are taken into account. However, 
linear-independence of all the functions with unknown parameters on the driving signal is 
crucial to a successful parameter identification. 

6 Conclusion 

In summary, concrete examples showing possible occurrence of failed parameter identifi- 
cation have been numerically given in the paper. The mechanism inducing this failure 
has been further rigorously interpreted. It has been pointed out that chaotic property 
of driving system is not always crucial to achievement of parameter identification either 
in a mathematical argument or in a numerical experiment. Actually, it is not the chaos 
but the hypothesis [LIM] that guarantees a successful parameter identification based on 
adaptive synchronization techniques. However, making good use of chaotic property might 
easily lead to validity of hypothesis [LIM]. Apart from linear-dependence of functions, ap- 
proximate linear-dependence of functions with parameters pending for identification on 
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the synchronized orbit should be always avoided in numerical simulation and even in real 
application. 

Furthermore, complete synchronization via new adaptive coupling techniques in a class 
of polynomial systems where nonlinearity is not globally Lipschitz has been theoretically 
investigated by virtue of the LaSalle invariance principle. Also it has been rigorously veri- 
fied that every trajectory generated by the coupled systems is bounded. By these derived 
theoretical results, our newly proposed technique is convinced to be a rigorous and feasible 
approach for realization of complete synchronization and parameter identification in the 
Lorenz-like systems. Besides, adaptive coupling techniques are also imported to realize 
parameter identification in systems with time-delay. Those discussion also shows the great 
importance of the condition that functions with parameters pending for identification on 
the synchronized orbit should be linearly independent. 
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Fig. 1: (Color online) A successful complete synchronization between the Lorenz system (U) and 
^ by means of the adaptive design coupling. Here, system ^ possesses an asymptotically stable 
equilibrium E ~ (6.8313,6.8313,1.6667)"'" instead of the strange attractor. The variation of the 
driving signal with the response state are shown in (a) and the evolution of response state in the 
phase plane are depicted in (b). Here, = 15, St = 2 and all the initial values are simply chosen 
as x° = ^ 10, = 1, e° = 1 {i = 1, 2, 3). 
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Fig. 2: The variation of the error between the parameters qi and pi with time initiating from to 
60 with step-size 0.01 {i = 1,2,3). In particular, qi fails to identify the accurate value of pi. All 
the parameters and initial values for coupling systems are the same as those given in Fig.l. 
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Fig. 4: (Color online) The variation of qi [i — 1, 2) with time initiating from to 15 with step-size 
0.01 when its initial value qig is differently taken [g = a, 6, c, d, e). Here, Vj — 2, 5j — 1 and all the 
initial values are taken as — — 10, e° = 1, j = 1, 2, 3. . 
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Fig. 5: (Color online) The strange attractor produced by system ([5]) are, respectively, plotted in 
the X1-X2-X3 plane (a) and in the time-state-X4 plane (b). The chaotic property is verified by the 
Lyapunov exponent portrait (c), where the largest Lyapunov exponent Ai is above zero. 
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Fig. 6: The variation of the error between the parameters qj and pj with time initiating from to 
25 with step-size 0.01. Indeed, 93^4 fails to identify the accurate value of ^3^4, respectively. Here, 



rj — 15, Sj ~ 2, and all the initial values are taken as Xj 



1, y°i = 6, 2/2° 2/3° = 10, 2/0 = 1, 



90 = 0, and e°^l (j = 1,2,3,4). 




Fig. 7: (Color online) The variation of X3 and 0:3 (1 
time, respectively 



on the synchronized orbit x*{t) with 
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Fig. 8: (a) Successful complete synchronization and parameter identification for chaotic driving 
signal when a = —2 and 6 = 4; (b) Failed parameter identification when a = 2 and 6=1. This 
failure is simply due to an approximate dependence between f(^x(t)^ and g(^x{t — t)^ in a macro 
scale. Here, both / and g are taken as sinusoid functions, time-delays are taken as t = 10, 6 — 2, 
and time step size is 0.01. 
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